Local current injection into mesoscopic superconductors 
for the manipulation of quantum states 



M. V. Milosevic, 1 A. Kanda, 2 S. Hatsumi, 2 F. M. Peeters, 1 and Y. Ootuka 2 

1 Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
2 Institute of Physics and TIMS, University of Tsukuba, Tennodai, Tsukuba, Ibaraki 305-8571, Japan 

(Dated: November 6, 2009) 

We perform strategic current injection in a small mesoscopic superconductor and control the 
(non) equilibrium quantum states in an applied homogeneous magnetic field. In doing so, we realize 
a current-driven splitting of multi- quanta vortices, current- induced transitions between states with 
different angular momenta, and current-controlled switching between otherwise degenerate quantum 
states. These fundamental phenomena form the basis for discussed electronic and logic applications, 
and are confirmed in both theoretical simulations and multiple-small-tunnel-junction transport mea- 
surements. 

PACS numbers: 73.23.-b, 74.78.Na, 85.25.Hv 



Preparation and manipulation of discrete quantum 
states are crucial for applications of quantum physics in 
nanoscale electronics, particularly switching devices and 
memories. The needed quantum states with discrete con- 
ductance levels, and the external control of those, have 
been recently realized in solid electrolytes [1], graphene 
2[ (both voltage biased), light-driven molecular switches 
3(], and magnetic nanowires tuned by magnetic field [4|. 
With a similar goal, the current-driven processes started 
to attract immense attention since the demonstration 
that the magnetization state of a nanomagnet can be 
influenced directly by electrical current [5j . The key con- 
venience is that locally applied current enables switching 
of individual submicron elements even in integrated elec- 
tronic circuits. 

The latest example of the fundamental role of cur- 
rent injection in low-dimensional physics is the achieved 
electric flipping of the magnetic vortex core [6], with- 
out change in the vortex chirality. Vortices, rotational 
flow of currents or matter with a characteristic cavity 
at the center, can be found in many subfields of quan- 
tum physics. In small magnetic elements, vortices owe 
their existence to shape anisotropy and demagnetizing 
fields 0- In semiconducting quantum dots the rotation 
of electrons is induced by external magnetic field and 
vortices may form if this rotation is strong [8]. Sponta- 
neous vortex nucleation has been recently observed even 
in exciton-polariton condensates [9[, a composite boson 
system resembling superfluidity in Bose-Einstein conden- 
sates in which vortices are readily found [lQj |. Still, vortex 
matter is most intrinsic to superconductors in a magnetic 
field. It is particularly rich in samples comparable to 
the coherence length £ and/or penetration depth A [Tlj . 
where vortex states become strongly influenced by sam- 
ple geometry [12]. Due to strong lateral confinement, 
vortices may even merge into a multi-quanta ("giant") 
vortex, otherwise unstable in an open geometry [13|, 
Such vortices contain multiple phase change of 2tt en- 
circling the central cavity, and are also found in rotated 



superfluids and in strongly confined quantum dots in a 
magnetic field [151 ]. 

Vortices not only represent a key feature of topologi- 
cally confined quantum-mechanical systems but may also 
be of use for electronic device applications. It has been 
shown recently that each change in the core of magnetic 
vortices gives rise to a measurable electric signal [l6| . 
Changes in vortex configurations in quantum dots intu- 
itively lead to distinct features in electronic transport 
[l?]]- In superconductors, the quasiparticle excitations 
inside vortices [18[ form coherent quantum-mechanical 
states that depend on the number and arrangement of 
trapped flux quanta. As a consequence, the sample con- 
ductance measured along the vortex lines is determined 
by the transparency of the given vortex configuration 
[191 ]. Therefore, with changing magnetic field and thanks 
to the consequent change in vortex states, different levels 
of conductance can be achieved and a quantum switch 
is realized. Up to date, varied magnetic field remains 
the only tool for control of the angular momentum of 
vortex states and their configuration. In this Letter, we 
present a new method for vortex manipulation in small 
superconductors - by strategically applied current. We 
demonstrate the first electronically triggered splitting of 
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FIG. 1: SEM image of the sample - the mesoscopic super- 
conducting square in a transport structure, with two tunnel 
junctions on top, as illustrated on the right (indicated direc- 
tions of field and current are denoted positive). The details 
of fabrication can be found in Ref. [141 ] . 
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the giant vortex as a fundamental novelty, controlled 
transitions between vortex states of different angular mo- 
mentum for a given magnetic field (i.e. a current-driven 
quantum switch) , as well as the manipulation of degener- 
ate vortex equilibria from one to another, useful for logic 
applications. 

Our sample is a d = 40 nm thick, 1.1 x 1.1 jim 2 Al 
superconducting square designed for the multiple-small- 
tunnel-junction (MSTJ) measurement [141 ]. in which two 
normal- metal (Cu) leads (width 0.25 fim) are connected 
to the top side of the sample through tunnel junctions 
(see Fig. [I]). In addition, the square has three direct Al 
leads (width 0.3 fim) centered at remaining sample sides, 
two of which are used for current (/) injection. The co- 
herence length was estimated to £(0) = 150 — 190 nm, and 
the superconducting transition temperature T c was 1.35 
K. The magnetic field sweep rate was 3 mT/min. In the 
MSTJ measurement, a small constant current (typically 
1 nA) is applied to both junctions and the junction volt- 
ages, V\ and V2, are measured simultaneously as a func- 
tion of the magnetic field or the injected current /. An 
observed voltage reflects the local density of states and 
also the local super current density underneath the junc- 
tion, because the energy gap decreases with increasing su- 
percurrent density. As a result, voltage jump corresponds 
to a transition between different vortex states. Further- 
more, the comparison of voltage signals at two junctions 
indicates the symmetry of the circulating current, and 
thereby V\ 7^ V2 points out the multi-vortex configura- 
tion, while V\ = V2 suggests a giant-vortex state. 

To characterize the stationary and dynamic properties 
of the device, we solved numerically the time-dependent 
Ginzburg-Landau (TDGL) equation [2(j 

7rrlw^4^ = (v-*a)V +( i-h> 

(i) 

coupled with the equation for the electrostatic potential 
Atp = div (3 : ('0*(V — iA)^)). Here, the distance is mea- 
sured in units of the coherence length £, ip is scaled by 
its value in the absence of magnetic field ^0, time by 
tql = 2Tft/7r0o, vector potential A by c?i/2e£, and the 
electrostatic potential by cpo = Ti/2erGL- T = 2t#V ; o/^ 
with te being the inelastic electron-collision time (for Al 
samples, te ~ 10 ns gives Y « 10 3 ). Parameter u = 5.79 
is taken from Ref. [20[. Note that in Eq. (pQ) the screen- 
ing of the magnetic field is neglected, as our sample is suf- 
ficiently thin to exhibit effective type-II behavior (there- 
fore A = \B{y^ — x, 0)). The current leads (see Fig. [1]) 
were simulated by imposing — V<£ = ji, where ji is the 
injected current density in units of jo = c<I>o/87r 2 A 2 £, 
A = X 2 /d. At the remainder of the sample edges, Neu- 
mann boundary condition was used (j± = 0). 

Fig. [2] shows the measured voltage and calculated 
Helmholtz free energy of the sample as a function of 
applied magnetic field {T — y j v \^\ A dV ^ in units of 
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FIG. 2: Measured voltage (a) and calculated free energy (b) 
at T = 0.7 K. Insets show phase contourplots of the supercon- 
ducting order parameter for L = 3 states (multi- and giant 
vortex) marked by open dots in free energy curves. 



JF = i^ 2 /87r, where V is the sample volume). In Fig. 
[2Kb) different energy levels correspond to states with dif- 
ferent angular momentum, i.e. vorticity L. Here taken 
£(0) = 152 nm and T c = 1.34 K provided the best agree- 
ment with experimental vortex penetration/expulsion 
fields [see Fig. [2]. As already pointed out, some vor- 
tex states may exhibit ring-like configuration of individ- 
ual vortices [13[, i.e. a multi- vortex [dashed lines in Fig. 
[2](b)], while in others all flux quanta coalesce into a sin- 
gle giant- vortex (solid lines). At higher temperatures, the 
characteristic lengthscales grow and effective confinement 
favors giant- vortex states. 

The distinction between the giant- and multi-vortex 
states in mesoscopic superconductors can also be made 
by monitoring of the vortex expulsion field B e for a given 
state as a function of temperature, as detailed in Ref. [2l[ 
(decreasing B e (T) tendency is found for a multi- vortex 
vs. increasing B e (T) behavior for a giant- vortex state). 
Our measurements and calculations indicated that for 
e.g. L = 3 vortex state, the giant-vortex form is fa- 
vorable above 0.9 K. In what follows, we investigate the 
influence of applied current on the shape of the vortex 
state. Following the scheme shown in Fig. [TJ we applied 
weak dc current in the direction down-to-right, such that 
vortices experience Lorentzian-type of force (F = x 3>) 
approximately towards top-left corner of the sample. Due 
to the spatial inhomogeneity of the force and the interac- 
tion with the sample boundaries, it can be energetically 
favorable for the giant-vortex to separate into individ- 
ual ones and accommodate the imposed conditions. We 
achieved this giant-to-multi vortex splitting in both the- 
ory and experiment at temperatures above 0.95 K and 
this is the first such result in vortex physics. As shown 
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FIG. 3: Experimentally and theoretically obtained vortex ex- 
pulsion field as a function of temperature for the L — 3 vortex 
state without and with applied current. Different behavior of 
the curves indicates multi- or giant- (shaded areas in experi- 
mental curves) vortex state. 



in Fig. O for e.g. T = 1.0 K, we found a giant- vortex 
in the absence of applied current, but multi- vortex when 
even a weak current of 5 jaA was applied. In supplemen- 
tary material [23], we also show the simulated dynam- 
ics of the giant-vortex splitting under applied drive at 
T = 1 K. Note however that in the case of opposite cur- 
rent (J = — 5 /iA), the drive 'drags' vortices towards the 
bottom-right corner, where superconductivity is strong. 
This compresses the previously existing multivortex, and 
favors the giant-vortex state (see corresponding experi- 
mental curve in Fig. [3] for 0.7 < T < 0.95 K). In ad- 
dition, the vortex expulsion field is higher compared to 
1 = case, since the specific Lorentzian drive restricts 
access of vortices to the sides of the sample. 

The above splitting of a giant-vortex by injected cur- 
rent is a non-equilibrium process, as giant- and multi- 
vortex realization of the same L-state cannot coexist. 
However, for a given multi- vortex state, there can be 
maximum n energetically degenerate vortex configura- 
tions with respect to n sample boundaries. This fact 
forms a base for the fluxonic cellular automata [22[ , where 
e.g. two vortices in a square can occupy two alternative 
diagonal positions, forming states that can be labeled as 
logic '0' and '1'. One of the essential problems in cellu- 
lar automata is the successful preparation of inputs, i.e. 
the feasible manipulation of bits, and this can be real- 
ized by local current injection. In Fig. 0J we show the 
manipulation of the L = 3 multi- vortex state at T = 0.7 
K. For sufficiently strong negative current, two vortices 
are attracted to the bottom-right corner rather than just 
one, which rotates the vortex structure by c.a. 30°. In 
other words, the applied current enables the transition 
between otherwise degenerate states, and the switching 
is fully reversible. This is demonstrated in Fig. [4j where 
we realized multiple switching back and forth between 
two degenerate states by an alternating current. The 
rotation of the vortex configuration is verified by the al- 
ternating voltages measured at the two tunnel junctions. 
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FIG. 4: Switching between degenerate vortex states (see in- 
sets) by alternating current (|/| <> —7.5 /jlA), demonstrated 
experimentally through measured voltages at junctions 1 and 
2 (see Fig. [Q). 



Obviously, here reported phenomena are sensitive to 
the sample parameters, but crucially depend on the mag- 
nitude of the applied current. The larger current exerts 
larger force on vortices, which may expel some of them 
and/or help nucleate new ones in cases when, roughly 
speaking, local current on one of the sample edges ex- 
ceeds the depairing current. The first scenario is partic- 
ularly interesting at applied magnetic field which stabi- 
lizes several vortex states (see e.g. Fig. [2j for B = 12 
mT). In such a case, one can induce controlled transitions 
between states of different vorticity by applying current 
of needed magnitude! We show this in Fig. Ufa), where 
one vortex is pulled out of the ground-state L = 3 giant- 
vortex (by / < 0), and double vortex [L = 2) remains 
stable after the current is switched off. Successful switch- 
ing manifests as a peak of voltage V vs. time. If the 
same state was exposed to even larger current, vortices 
are sequentially expelled from the L = 3 state, but also 
new vortices appear at the opposite edge of the sam- 
ple. Eventually, a 'dynamic equilibrium' is established 
(as shown in Fig. 03(b) and in supplementary animations 
[231]). exhibiting clear periodicity in the V(t) signal linked 
to entry and exit of vortices. Latter oscillations of volt- 
age are tunable by applied current, but in the present 
(Al) sample they are mostly in the GHz frequency range 
{tgl « 0.2/(1 - T/T c ) ps). 

In Fig. [5](c) we show the experimental evidence for 
current-driven transitions between states of different an- 
gular momentum, at T = 0.85 K [24]. After decreasing 
magnetic field to B = 11.2 mT, the L = 4 to L = 3 tran- 
sition was induced by applied current of / = — 15 /iA, 
analogously to the scenario shown in Fig. 03(a). Voltage 
measured in further decreasing field shows that L = 3 
state indeed remains stable. Therefore, our sample can 
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FIG. 5: Calculated voltage between current leads vs. time 
(at T — 1.0 K), for the case of (a) pullout of a flux-quantum 
(transition from L — 3 to L — 2) , and (b) active vortex pene- 
tration and expulsion, with average vorticity L — 3 (insets de- 
pict snapshots of the Cooper-pair density, see corresponding 
animations in [23]). (c) Experimentally measured voltage vs. 
decreasing applied field, prior and after the current-induced 
L — 4 — > 3 transition. 



act as an excellent quantum switch, where applied cur- 
rent opens/closes the ballistic channel for quasiparticles 
through the vortex core (see Ref. 0). 

Latter L = n <-> n + 1 transitions can be achieved 
for both polarities of the applied current, as one fluxon 
can be either 'pushed' or 'pulled' out of the sample. In 
Fig. [6] we give the experimentally obtained values for 
the threshold current needed for enforced transitions be- 
tween states with vorticity 0<L<4atT = 0.85K, both 
for increasing and decreasing magnetic field. The typical 
values of minimal current needed for the switching were 
in the |/| < 30 fiA range for T = 0.85 K, and grew sig- 
nificantly with decreasing temperature (e.g. |/| < 80 fiA 
range for T — 0.7 K). Note that in some cases it becomes 
irrelevant if we perform vortex penetration or expulsion 
- e.g. at B = 7 mT, applied current of / > 30 fiA will 
provide L = 1 state, regardless of the initial vortex state 
(L = or L = 2). 

To summarize, we realized electronic control of the 
quantum states in a mesoscopic superconductor by 
strategic injection of dc current. Out of equilibrium, we 
performed the splitting of a giant, multi-quanta vortex 
into individual vortices. In equilibrium, we demonstrated 
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FIG. 6: Minimal switching current for states with < L < 4 
(experiment done at 0.85 K), with distinction between the 
L — > L — 1 and L — > L + 1 processes (maximal current in 
experiment is limited by the critical current of the leads). 
The quantum switch was not operational in the shaded areas. 



the current-induced transitions between states with dif- 
ferent angular momenta, which also has a step-like effect 
on local ballistic conductance and is applicable as a quan- 
tum switch. Finally, switching between dual equilibria 
was also achieved, most useful for the fluxonic cellular 
automata. Local current injection therefore proves to be 
a powerful method for manipulation and preparation of 
quantum states on a submicron scale, and further work 
is needed to test the analogies in other quantum systems. 
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